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The spin correlations w^, r = 1, 2, 3, and the probability pn of finding a system in the Neel state 
for the antiferromagnetic ring Feg 11 (the so-called 'small ferric wheel') are calculated. States with 
magnetization M = 0, total spin < S < 15 and labeled by two (out of four) one-dimensional 
irreducible representations (irreps) of the point symmetry group Dg are taken into account. This 
choice follows from importance of these irreps in analyzing low-lying states in each 5-multiplet. 
Taking into account the Clebsch-Gordan coefficients for coupling total spins of sublattices (Sa = 
Sb = t^) the global Neel probability pjj can be determined. Dependencies of these quantities on 
state energy (per bond and in the units of exchange integral J) and the total spin S are analyzed. 
Providing we have determined Pn(S) etc. for other antiferromagnetic rings (Feio, for instance) we 
could try to approximate results for the largest synthesized ferric wheel Feis- Since thermodynamic 
properties of Fe6 have been investigated recently, in the present considerations they are not discussed, 
but only used to verify obtained values of eigenenergies. Numerical results are calculated with high 
precision using two main tools: (i) thorough analysis of symmetry properties including methods of 
algebraic combinatorics and (ii) multiple precision arithmetic library GMP. The system considered 
yields more than 45 thousands basic states (the so-called Ising configurations), but application of 
the method proposed reduces this problem to 20-dimensional eigenproblem for the ground state 
(5* = 0). The largest eigenproblem has to be solved for S = 4; its dimension is 60. These two 
facts (high precision and small resultant eigenproblems) confirm efficiency and usefulness of such an 
approach, so it is briefly discussed here. 
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I. INTRODUCTION 

Magnetic hysteresis with a pure molecular origin first 
time was observedEltl in the mixed-valence Mn complex 
[Mn 12 Oi 2 (CH3COOi 6 (H20)4]-2CH 3 COOH-4H 2 (com- 
monly referred to as Mni 2 .acetate), synthesized and in- 
vestigated by Lis in 1980.B Since this discovery meso- 
scopic magnetic systems, comprising transition-metal 
ions, have attracted a great deal of attention. Many 
large molecular aggregates have been fabricated and in- 
vestigated both experimentally and theoretically (includ- 
ing numexkjal simulations) . This includes. Mn molecules 
like MngJj'tJ a ferromagnetic copper ringjj mixed Mn-Cr 
compoundsLhke [Cr ln (CNMn n L) 6 ](C104)9,Q some nickel 
complexes, au and many polynuclear iron clusters-known 
as Fe n with n varying from twaJ through 20.Ej'Ell Some 
interesting features have been also observed in, a crystal 
made of nanometric molecules (so-called Vi5)E3 and in a 
single cobalt nanocluster .Emilia 

Such systems of mesoscopic dimensions containing 
transition-metal ions and oxygen are attracting increas- 
ing interest for several reasons, including the possibil- 
ity of observing quantum phenomena on a macroscopic 
scale. It may lead to better understanding of the mecha- 
nism associated with the transition from simple paramag- 
netic to bulk magnetism (also low-dimensional in mag- 
netic chains). On the other hand, clusters comprising 
transition elemeats exist in several metalloenzymeSpaBd. 
metalloproteinala as ferritin storing Fe in mammals. t3'tZI 



Moreover, nanostructures are part of the_development of 
high densite, magnetic recording mediaEfl and quantum 
computing.!!^ 

Molecular magnets show both antiferro- and ferromag- 
netic behavior, though in m&ny cases single spins are cou- 
pled antiferromagneticallylij Incomplete compensation 
of this coupling may lead to the large s pin gr o u nd, state of 
a cluster, e.g. S = 10 in Fes and Mni 2 ,E3'li3i2£lO'E2l what is 
related to a strong easy-axis term (Ising- type anisotropy). 
There are also purely ferromagnetic molecules as dodc- 
canuclear nickel complexJjNi^) with the ground state 
characterized by S — 12.BErtL3 The other class of nano- 
magnets comprises molecules where the cancellation is 
complete and the ground state has S — 0. The most 
profound example is provided by a series of antiferro- 
magnetic Fe(III) rings Fe@, Feio, Fei 2 , and Fei 8 £j Re- 
cently a new octanuclear-.cluster Fes has been synthe- 
sized and investigated rW^\ A deviation from coplanarity 
in the case p£Eeio ring, about 9.8 A in diameter, is less 
than 0.01 AJlj'Ej so these rings can be considered as two- 
dimensional structures of very high symmetry: at least 
idealized or approximate point-group symmetry isomcp=. 
phic to the dihedral group D„ (Sq or Ds^, for example).E£l 

Due to the presence of organic ligands wrapping the 
clusters, the inter-cluster interactions are vanishingly 
small (or can be reduced] by dissolving the clusters in 
the appropriate solvent). t2l Each magnetic object has per- 
fectly defined size and new techniques of molecular ch om -. 
istry provide ensembles of iso-oriented nanomagnets.E2lEZl 
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Properties of such systems can be very well simulated 
using contemporary-,|Gomputers within the framework of 
finite spin modelsflEa However, direct calculations of 
eigenenergies and eigenstates, even using the Lanczos ex- 
act diagonalization technique, are subjected to the so- 
called "combinatorial explosion" . For example, in the 
case of the largest "ferric wheel" Fei 8 one has to deal with 
more than 10 14 basic states. This is especially important 
in the case of antiferromagnetic systems when the ground 
state lies in the subspace of states with the total magneti- 
zation M = and determination of some low-lying states 
only is a formidable task. This can be done a bit more 
manageable exploiting symmetry properties of a given 
Hamiltonian 7i (including the time-reversal symmetry) 
and good quantum numbers, i.e. eigenvalues of operators 
commuting with TL. In the latter case we consider the 
total magnetization M (the eigenvalue of S z = Y^j=i s j) 
and the total spin S (the eigenvalues of the squared total 
spin S 2 = (Y^+i s j) 2 are S(S + 1), of course). While 
model Hamiltonians commute with S z as a rule, the acL 
ditional anisotropy terms break commutation with S 2 .EZI 
However, if the isotropic Heisenberg interactions 



(1) 



(the sum is taken over the nearest-neighbor pairs) are 
dominant, then the effect of magnetic anisotropy and 
Zeeman interactions, can be considered separately for 
each 5-multiplet.EZrEa Therefore, in this article we restrict 
ourselves to Hamiltonians commuting with S 2 . It should 
be stressed that this includes interactions with different 
ranges (e.g. next nearest-neighbors) and forms (e.g. bi- 
quadratic terms (Sj ■ s,t) 2 ). 

As regards symmetry properties, the simplest task is 
to take into account the time-reversal symmetry. It is 
enough to consider states with M > since the oth- 
ers (with M < 0) are determined in a straightforward 
way. The largest eigenproblem dimension is reached con- 
sidering the subspace L comprising states with M = 
(the are more than 5.5 • 10 12 such basic states for s = | 
and n = 18). The point-group symmetry G enters the 
problem considered via the Schur Lemma: eigenspaces 
of operators commuting with all representation opera- 
tors P(g), g G G, are labeled by the irreducible repre- 
sentations (irreps) of G. However, numbers n(T) in the 
decomposition 



p = 0n(r)r 



(2) 



are very large, even if P is restricted to Lq. Let L Vr , 
1 < r < n(T), be copies of the eigenspace of an irrep T 
spanned over vectors {^7) | 1 < 7 < dimi r }. Solv- 
ing an eigenproblem for any operator of symmetry G we 
have to construct an n(T) x n(r) matrix for an arbitrarily 
chosen, but fixed, vector index 7. To do it in an efficient 
way, the single, meaningless, repetition index r has to be 



replaced by a series of indices with well-determined math- 
ematical notion. In the case considered the represen ta- 
tion P is in fact a permutation representationjEaE^lO so 
the methods of algebraic combinatorics give an efficient 
way to deal with this problem. Such approach to finite 
spin systems has been discussed recently in a series of 
papersE3'E3 In this work we aim to present and discuss 
results obtained for the "small ferric wheel" Fe§. The 



reason of such a choice is 
has been well investigated,! 



(i) this molecule 
so we can verify 



our results; (ii) the number of states are relatively small 
(6 6 = 46 656 and dimio = 4332), then many character- 
istics can be obtained without tedious calculations; (iii) 
all numerical algorithms and procedures can be tested 
before going to larger problems. 

The paper is organized as follows. In Sec. [n] a short 
account of the iron antiferromagnetic rings is presented. 
Combinatorial classification of finite spin system states 
is briefly discussed in Sec. |l|. Section IV contains 
eigenproblem solutions for Fe 6 cluster. Eigenstates ob- 
tained are analyzed in Sec. [v| (e.g. static spin correla- 
tions are calculated) . Some concluding remarks are made 
in Sec. |Vl| Due to reasons given in Sec. IV states cor- 



responding to two one-dimensional irreps are considered 
only. Hence, the energy spectrum is not determined and 
the thermodynamic quantities are not calculated. 



II. ANTIFERROMAGNETIC FERRIC WHEELS 

There are several mesoscopic clusters comprising iron 
ions bridged by oxo-groups with different magnetic prop- 
erties. Some of them have the ground state with high 
total spin S, e.g. Fei II (OCH 3 ) 6 (dpm) 6 j&jith S = 5,1^1 
[Fe I 8 II (tacn)0 2 (OH) 12 ]Br 8 rwith S = lojil^ and ferro- 
magnetic Fe^ 1 with S — 8.E3 On the other hand, in some 
systems antiferromagnetic coupling results in complete 
cancellation of local spins leading to nonmagnetic (S = 0) 
ground state; all.such molecules contain iron(III) ions a: 
ranged in cyclesEiJ or dimers as in [Fe(OMe)(dbm)2]2 
In the first case we have three types of Fe6 riugs, a 
medium ferric wheel [CsCFe 8 {N(ClLGH20)3}8]ClJ 2 l the 
Feio(OCH)3)2o(C*H 2 2 Cl)io ringM and the largest 
ferric wheel Feig.E2l There is also a ring-shaped molecule 
Fei2 with a net spin S = in the gepund state, but it does 
not exhibit full cyclic symmetryJlJ The number of basic 
states is given as 6™, n = 2,6,8,10,12,18, so we deal 
with eigenproblems of dimensions 36, 46 656, 1679 616, 
60 466176, 2176 782 336, and 101559 956 668 416, respec- 
tively. The dimensions of Lq, the subspace containing 
states with M = are as follows: 6, 4332, 135 954, 
4395456, 144 840476, and 5542 414273 884; the approx- 
imate ratios 6"/dimL are 6.00, 10.77, 12.35, 13.76, 
15.03, and 18.32, respectively, so dimLo increases a bit 
slower than 6". 

In this work the small ferric wheel is considered; the 
following thtee molecules of this type are known and 
investigated:^ 
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[NaFe 6 (OCH 3 )i 2 (C 17 H 15 04)6]C104 , 
[NaFe 6 (OCH 3 )i 2 (Ci5H 11 02)6]C104 , 
[LiFc 6 (OCH 3 )i 2 (C 1 5H 11 2 ) 6 ]C104 . 

The first cluster has crystallographically imposed Sq 
point-group symmetry, whereas in the other cases this 
symmetry is approximated Therefore, the results ob- 
tained here are mostly relevant to the first one. 

III. CLASSIFICATION OF STATES BY MEANS 
OF ALGEBRAIC COMBINATORICS 

The method proposed can be roughly divided into four 
steps: combinatorial, linear, magnetic, and numerical. 
The last one consists in efficient implementation of pro- 
posed algorithms, working-out methods of storing of con- 
figuration etc. These, problems have been recently pre- 
sented elsewhere,Eiloc3 so they are left out in the present 
paper. The second and third steps are slightly modified 
standard procedures known since the famous paper of 
Bonner and Fishercil and reused lately in many articles, 
for example by Waldmann in Ref. Therefore, below 
we concentrate on the first step since combinatorial meth- 
ods (especially methods of algebraic combinatorics) are 
not so popular. More detailed discussion on this topic, 
from the mathematicaL-point of view, can be found in a 
monograph by KerberO 

To begin with we have to introduce the notion of a 
(finite) group action and define two structures: an orbit 
and a stabilizer. We say that a group G acts on a set X if 
for each g £ G and x £ X a product gx £ X is uniquely 
determined and for all g, g' £ G, x £ X 

(gg')x = g(g'x) and ecx = x , 

where &q denotes the unit element in G. If both G and 
X are finite we call the action finite. An orbit G(x) C X 
of an element x £ X is a set containing all x' £ X which 
can be obtained from x, i.e. 

G(x) = {gx \g£G}. 

Since an orbit is an equivalence class then it can be rep- 
resented by any of its elements. For each x £ X its 
stabilizer 

G x = {g | gx = x} C G 

contaips elements of G leaving x invariant. It can be 
shownEa that orders |G(x)| and \G X \ of an orbit and its 
stabilizer, respectively, are related by the following for- 
mula 

\G\ = \G(x)\\G x \. 

Moreover stabilizers of all elements in a given orbit G(x) 
are conjugated to each other (as subgroups of G); a class 
of conjugated subgroups determines a type of an orbit. 



For any set Y the action of G on X can be r.ai3ed to an 
action of G on functions /: X — > Y as followsEa 

Vx£X {gf){x) = f{g- 1 x). (3) 

In the case of finite systems the so-called Ising con- 
figurations (basic states) can be considered as mappings 
fi: X — ► Y , where X = {1, 2, . . . , n} is a set of node la- 
bels and Y — {— s, — s + 1, . . . , s} is a set of ^-projections 
for the spin number s. Basic states |mi, to 2 , . . . , m n ), 
rrij £ Y for 1 < j < n, are in a one-to-one correspon- 
dence with functions (i providing that for all j £ X the 
equality fj,(j) = nij holds. The largest (discrete) group 
acting on X is the symmetric group £„ comprising all n\ 
permutations a of the n-element set X. According with 
Eq. (||) one obtains 

{<J l x){j)= i x{a- 1 j) (Vl<j<n) 

or, writing it for basic states \m±, m 2 , . . . , m n ), 

a\mi,m 2 , . . . ,m„) = \m <7 -i 1 ,m a -i 2 , ■ ■ ■ ,m a -i n ) . 

Let us consider, for example, a system of n — 4 spins 
s = |. One of the Ising configurations is a state \ip) — 
|o> 2> |j §)• Acting with all 4! = 24 permutations on this 
state one obtains only six different states (i.e. |G(|i/'))| = 
6), since a stabilizer Gy, is given as 

G* = {1, (12), (34), (12) (34)}, 

where (jk) denotes a transposition of nodes j, k and 1 is 
the unit element (identity) in E ra . To obtain all elements 
of the orbit G(\ip)) it is enough to act with representa- 
tives of left cosets Y^^jG^,, for example with permutations 
1, (13), (14), (23), (24), and (14)(23). If a is a transposi- 
tion then cr _1 — a, so 

G(UIM — /|I 1 3 3\ |3 1 1 3\ |3 1 3 1\ 
"Al *r I ) U2'2'2'2''12'2'2'2''l2'2'2'2'' 
|I 1 I 3\ |1 3 3 1\ |3 3 1 l\i 
l 2 > 2' 2' 2 /'l 2 ' 2 ' 2 ' 2 /'l 2 ' 2 ' 2' 2' > ' 

Note, that this orbit is determined in unambiguous way 
by a nonordered partition [0, 0, 2, 2] of n — 4 into 2(3/2) + 
1 = 4 non-negative parts. Each entry ki, I = 0, 1, . . . , 2s, 
of [k] = [ko, fci, . . . , /c 2s ] denotes a number of projections 
Til j , j = 1, 2, . . . , n, equal to I — s. Since [k] is a partition 
of n, then ki = n and fc; > 0. An orbit containing 

states with ki projections I — s is denoted hereafter as 
0[k] and is represented, for example, by a configuration 

I ~s, s,. . . , -s , -s + 1, . -s + 1 , . . . , s,s,...,s) . 

ko times k\ times k-2 s times 

Moreover, a nonordered partition [k] determines the mag- 
netization M for all states in 0[k] since 

2s 2s 

M = - s)ki = ^2lki-ns. (4) 

1=0 1=0 



4 



In the example presented one obtains M = — ^-O— §-0+|- 
2 + | • 2 = 4. The same magnetization M = 4 is obtained 
for states from the four-element orbit O[0, 1,0,3] repre- 
sented by a state | — \, §, §, §). Some group-theoretical 
properties of orbits 0[k] depend on their type, i.e. on 
a class of conjugated subgroups of S n a stabilizer of 
states in 0[k] belongs to. Such classes are reocesented 
by the so-called Young subgroups £r re i C £„,ESEJ where 
[k] is an ordered partition [kq, k\, . . . , k z ], z < 2s + 1, 
of n into no more than 2s + 1 nonzero parts ko > 
Ki > . . . > k z > 0. The Young subgroup £r re i is a di- 
rect product of symmetric groups £ Ki , I = 0,1,..., z. 
Each group £ K( contains permutations of a K;-element set 
{Ki + l,Ki + 2,...,K + Kl }, where K x = «i- The 

nonordered partition [0,0,2,2] is of a type represented 
by the ordered partition [2, 2] with a number of nonzero 
parts equal to z + 1 — 2. Therefore, a class of conjugated 
orbits is represented by the Young subgroup £2 <8> £2, 
where the first group £2 contains permutations of the set 
{1,2} (K = 0) and the second factor £ 2 — of the set 
{3, 4} (Ki = 2). Other orbits of the same type are repre- 
sented by states \a, a, b, b), where a,b = — |, —5, i, | and 
a 7^ 6. A relation between nonordered partitions [k] of 
type [2, 2] and states representing orbits 0[k] is presented 
below; the magnetization M for states in a given orbit is 
also calculated: 



nonordered partition 


orbit 


representative 


M 


[2,2,0,0] 


1-3/2,- 


-3/2,-1/2,-1/2) 


-4 


[2,0,2,0] 


1-3/2,- 


-3/2, +1/2, +1/2) 


-2 


[2,0,0,2] 


-3/2,- 


-3/2, +3/2, +3/2) 





[0,2,2,0] 


l-l/V 


-1/2, +1/2, +1/2) 





[0,2,0,2] 


1-1/2,- 


-1/2, +3/2, +3/2) 


2 


[0,0,2,2] 


I + 1A- 


hl/2,+3/2,+3/2) 
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In the actual physical problems a symmetry group G 
is, or at least can be embedded subgroup of £„ . In a 
general case a restriction £„ J, G yields a decomposition 
of an orbit 0[k] into orbits of G. This decomposition de- 
pends on a type [n] of a nonordered partition, not on [k] 
itself. In the example presented each six-element orbit 
of type [2, 2] (determined by the action of £4) is decom- 
posed into two orbits of D4 C £4. The first is represented 
by the state \a, a, b, b) and contains four states, since its 
stabilizer (in D4) contains only the identity E and the 
two-fold rotation Up* which, as an element of £4, is a per- 
mutation (12)(34);Ej since the stabilizer is two-element 
subgroup of the eight-element group D4, then the orbit 
has 8/2 = 4 elements: \a, a, b, b), \a,b,b 7 a), \b,b,a,a), 
and \b, a, a, b). The other two elements, \a,b,a,b) and 
\b,a,b,a), form a two-element orbit with the stabilizer 
Dij = {E, G 2 , Uo, C/ 2 }. All six orbits presented above 
have the same decompositions into orbits of D4. More 
detailed discussion of combinatoric properties of the sys- 
tem of four spin s = I is presented in the Appendix. 

Therefore, we have determined a few starting steps of 
our procedure: (i) generate ordered partitions [k] of n 
into no more than 2s + 1 nonzero parts; (ii) find a decom- 
position of an orbit 0[k] into orbits of the Hamiltonian 



symmetry group G; (iii) for a chosen —ns < M <ns de- 
termine all nonordered partitions [k] satisfying the con- 
dition (|j); (iv) decompositions of orbits 0[k] into othita 
of G are analogous to those determined in Step (ii)Ed€3 
Note that orbits determined by the action of G are col- 
lected into types labeled by classes of conjugated sub- 
groups in G (in fact a representative U C G is used as 
such a label; this subgroup is a stabilizer of an element 
in an orbit under question). In this way an Ising state 
H can be labeled by the following indices: magnetization 
M, a partition [k], a stabilizer U, a representative v of 
an orbit G{v) 9 /1 (G v — U), and a representative (? r 
of a left coset g r U C G identifying /1 in the orbit G(v) 
{g r v — /i). These indices play different roles in the fur- 
ther considerations: if a given Hamiltonian TL commutes 
with S z then M is a good quantum number and can be 
used as an additional label of eigenspaces; U, v, and g r 
are used below to write down expressions for matrix ele- 
ments of spin operators (such as Ti. or S 2 ) in a compact 
and invariant form;EirE3 a partition [k] (together with an 
ordered partition [k]) is an additional index used during 
generation of Ising configuratio»s,|aHd it is very useful in 
calculations of matrix elements EUcIJ 

Now the linear structure of the space of states L comes 
into play. The permutation representation P restricted 
to an orbit of group G, and considered as a vector (lin- 
ear) representation of G, decomposes inta_irreps T in 
the same way for all orbits of a type [7.E3 For exam- 
ple, there are two orbits leading to M = and with 
D 2 C D4 being a stabilizer represented by ||, — |, |, — |) 
and |i, — 5, 5, — 5), respectively (n = 4, s = | as in 
the previous example). In both cases representation P 
restricted to an orbit is decomposed into a direct sum 
A\ ffi B\. For each orbit we construct the irreducible 
(symmetry adapted) basis 

\G\/\U\ 

|rVy>= £ aF r ^\g r v). 

r=l 

Since at this moment the considerations are limited to 
one orbit G(v), then indices distinguishing orbits (the 
magnetization M and the partition [k]) can be omitted. 
The formula presented is valid for all orbits G{v) of type 
U, i.e. for orbits represented by a state v with a stabi- 
lizer G„ = U. The index v distinguishes copies of an 
irrep T in the decomposition (^) restricted to an orbit 
G(v). Note that n(T) in such-sif-case is not larger than 
a dimension [r] of an irrep r.l2j't3 The coefficients a^. vl 
can be determined by the standard Haethods used in the 
case of permutation representations .Ej To obtain matrix 
elements of any operator H commuting with all P(g) we 
have to combine all vectors labeled by a given irrep T 
(including all copies). On the other hand, it follows from 
the Schur Lemma that it is enough to consider only one 
vector 7. It is rather tedious than difficult to derive a 
general formula for matrix elements Hxj vv ,u'v'v' labeled 
by orbit stabilizers U and U' , orbit representatives v^nd. 
v' , and indices v and v' distinguishing copies of foEl 
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The formula obtained contains products of two factors: a 
matrix element {Uvg r \H\U'v'ec) and a group-theoretical 
parameter of a model under discussion. Determination 
of the first one is a numerical problem, so it is not dis- 
cussed here— .More details can be found in other papers 
of authors .lEla The second factor depends only on the 
symmetry group G and can be determined once for all 
models with a given symmetry group. Moreover, it is 
possible to determine such factors in an analytical form 
for a family of groups, for example it has been done re- 
cently for the dihedral groups D„ describing symmetry 
of molecular rings £3 

The "magnetic" step of the procedure proposed can 
be applied to models with dominating Heisenberg term 
(H). At first we take into account the magnetization M 
restricting all considerations to subspaces Lm- We deter- 
mine matrices of the operator S 2 for all irreps T and solve 
eigenproblems S 2 |V>) = S(S + for < S < ns. In 

fact, since each S'-multiplet contains a state with M = 
(if ns is integer), then states \SM) can be obtained 
from states 1 5" 0) acting with the total step operators 
S ± = X}?=i ^ we are interested in the ground state 
of a bipartite antiferromagnet only, then we can limit 
calculations to the case S = and the one-dimensional 
irrep. In most cases this is the unit representation Tq 
(Ai in the case of dihedral groups). However, if s is half- 
integer and n is not divisible by four, then the ground 
state is labeled by the representation usually denoted 
as Ti (Bi for dihedral groups). This follows from the 
Marshall criterion:c3o Ising configurations related by 
the action of any pair of operators s'j's^ have coeffi- 
cients with opposite signs in the decomposition of the 
ground state into the basic states. In the case considered 
the ground state has to contain the Neel configuration 
|N1) = \s, — s, . . . , — s). According with the Marshall rule 
the other Neel state, |N2) = |— s, s, . . . , s), should enter 
the ground state with the same (opposite) sign if ns is 
even (odd). Hence, the ground state has to contain the 
state |Neel) = (|N1) ± |N2))/v / 2 which behaves as the 
basis vector of T {T\, respectively). 



TABLE I: Decompositions of orbits of the symmetric group 
£6, labeled by the partitions [k], into orbits of the dihedral 
group Dg, labeled by the stabilizers U. Each entry denotes 
number of orbits of type U contained in an orbit of type [n]. 



[k] 


[1,1,1,1,1,1] 


[2,2,1,1] [3,1,1,1] [3,2,1] [4,1,1] [3,3] 


[5,1] 


Ci 
D? 
Dg 


60 




14 10 4 2 1 
2 2 1 1 
1 



1 




[1,1,1,1,1,1], [2,2,1,1], [3,1,1,1], [3,2,1], [4,1,1], [3,3], and 
[5,1]; these (ordered) partitions represent types of orbits 
of the symmetric group SV There are 1, 14, 6, 4, 2, 3, and 
2 orbits (of states with M — 0) labeled by non-ordered 
partitions represented by each of the seven ordered parti- 
tions, respectively (together 32 orbits). For example, the 
ordered partition [3,1,1,1] represents six nonordered par- 
titions [fco, fci, . . . , ks] such that ^2 W = 15; these parti- 
tions and representatives of the corresponding orbit 0[k] 
are as follows: 



[1,0,1,3,1,0] 
[1,1,1,0,3,0] 
[0,1,3,1,0,1] 
[1,1,0,3,0,1] 
[1,0,3,0,1,1] 
[0,3,0,1,1,1] 



5 _1 1 1 1 3\ 

2 > 2 ' 2 ' 2 ' 2 > 2' ' 
5 _3 _ 1 3 3 3\ 
2' 2' 2 ' 2' 2 ' 21 ' 

3 _1 _1 _1 1 5\ 
2' 2' 2' 2 ' 2 ' 2' 
5 _3 1 1 1 5\ 

2' 2'2'2'2'2'' 
5 _1 _1 _1 3 5\ 
2' 2' 2' 2'2'2' 
3 _3 _3 1 3 5\ 
2' 2' 2 ' 2' 2 ' 21 ' 



Taking into account orbit cardinalities, determined as 
6!/fco!fci!fc2!^3!^4!fc5!, one obtains 

1*720+14*180+6*120+4*60+2*30+3*20+2*6 = 4332 

Ising configurations with M = what agrees with the 
number given above. 

We consider the nearest-neighbor interactions in a ring, 
so Eq. (0) can be rewritten as 



IV. EIGENPROBLEMS AND THEIR 
SOLUTIONS FOR THE SMALL FERRIC WHEEL 



3=1 



3 + 6 = j 



(5) 



The dimension of the subspace Lq in the case of six 
l ... i ■■>'->•) c;, i;,,. r — 4221 



spins s = | is equal to 4332. Since dimLi 



there are 111(= 4332 - 4221) states with S 
At first we determine such non-ordered partitions [k] 
[k , fci, . . . , k 5 ] of N = 6 into 2s + 1 = 6 parts that 

5 5 



1=0 



1 = 



There are 11 ordered partitions of N = 6 but only 
seven of them lead to non-ordered partitions satisfy- 
ing the condition presented above. An algorithm dis- 
cussed in Ref. 31 generates them in the following order: 



The— Hamiltonian point-symmetry group is G = Sg ~ 
D 6 .El There are 10 classes of conjugated subgroups for 
this group, but on% seven of them can play a role of 
an orbit stabilizer.E^I In the case considered, i.e. for the 
partitions presented above, only three subgroups appear 
in the decompositions of orbits of Eg: Ci, D^ and D§. 



These decompositions are presented in Table |. For ex- 
ample, each of six 120-element orbit of type [3,1,1,1] 
decomposes into ten 12-element orbits with the trivial 
stabilizer Ci. Taking into account numbers of orbits of 
Tjq one obtains: 339 orbits with a stabilizer Ci, 43 orbits 
with a stabilizer D°, and 3 orbits with a stabilizer Dg. 
These are all 385 orbits leading to the total magnetiza- 
tion M = 0. 



6 



TABLE II: Numbers n(T) for M = 0, 1 and S = 0; the last 
one is simply calculated as the difference of the previous two. 



r 




Ai 


A 2 


Si 


B 2 


Si 


E 2 


n(r), M 


= 


385 


339 


385 


339 


721 


721 


n{T), M 


= 1 


383 


325 


365 


334 


699 


708 


n(T), S 


= 


2 


14 


20 


5 


22 


13 



In the next step we take into account decompositions 
of transitive representations i? D6:C/ (i.e. the permutation 
representation _E restricted to an orbit of type U) into 
irreps T of D 6 .E3 The three stabilizers mentioned above 
give the following decompositions: 



RP 6 ' Ci = A 1 ®A 2 ®B 1 ®B 2 . 
RP e:T>1 = A 1 ® B 1 (B E 1 ® E 2 : 
jf>e:T>3 = A 1 ®B 1 . 



2E 1 2E 2 : 



This completes combinatorial and group-theoretical clas- 
sifications of states with M = 0. In a similar way one 
can classify states with M = 1, what enables us tq-deter- 
mine irreducible representations related to S = 0.E3 The 
numbers n(T) of subspaces with a given symmetry T for 
M = 0, 1 and S — are collected in Table 

Using the formulas discussed-, in Sec. [11], and pre- 
sented in the previous papers matrices of the Hamil- 
tonian H and the squared total spin S 2 can be con- 
structed. Of course, one can construct matrices for each 
r = A\,A 2 ,Bi,B 2 ,E\,E 2 . However, considering the 
ground state it is enough to consider T = B\. the ground 
state is a linear combination of 385 Ising states with 
M = transforming as |-Bi&i). We take into account 
also F = Ai sincfj. starting from a theorem proved by 
Lieb and Mattisj£3o it can be shown that these repre- 
sentations label alternately states with the lowest energy 
for S = 0, 1, 2, ... , 15 (At for odd and B x for even 5, 
respectively). The eigenstates determined for all 5* and 
r = A\,B\ are used to transform the Hamiltonian ma- 
trix to a block (quasi-diagonal) form. For example, to 
find the ground state and its energy we have to solve 
20-dimensional eigenproblem (cf. Table ||). 

A short note on numerical solutions are in place 
here. In the case of S 2 operator and one-dimensional 
irreps eigenproblems can be solved exactly. It has been 
achievedjusing the integer functions included in the GMP 
package.efl The square roots yielded by the action of 
operators _2Jid following from the formulas for ma- 
trix elements,E3 can be easily removed by the appropri- 
ate substitution. El In this way each eigenproblem for S 2 
can be written as a system of homogeneous linear equa- 
tions with integer coefficients. Therefore, the solution(s) 
can be expressed by elements of the smallest number 
field containing integers, so — in a general case — by ra- 
tional numbers also easily maintained with the use of 
the GMP package. The Hamiltonian matrix, with entries 
calculated with the GMP-functions, is transformed to a 



block-form using the float functions from this package.Ea 
The default precision is set to 256 bits. At this mo- 
ment a procedure for solving eigenproblems using the 
float GMP-functions has not been finished, so the eigen- 
values and eigenvectors are determined using the stan- 
dard procedures on 64-bit IRIX computers. The largest, 
60-dimensional, eigenproblem has been reached for 5 = 4 
and r = B\ . Obtained vectors are verified checking their 
scalar products (including norms) and their eigenvalues 
for S 2 and H matrices. 

In this way we have obtained 770 states expressed as 
linear combinations of the Ising configurations. Since 
states with 5 > represent in fact multiplets then ener- 
gies of 8436 states have been determined (4291 A\ and 
4145 B\ states). The ground state energy is equal to 
Eq = — 43.93471J; for other values of the total spin 5 
the lowest energies are in good agreement with Lande's 
rule A s = [E (S) - E<MJj= 5(5 + l)/3, which approxi- 
mates the energy gaps.EaEH The largest relative deviation 
(3.75%) is reached for 5=1 when Lande's rule gives 
Ai = |, whereas the value obtained is equal to 0.69169. 
The other gaps (5 = 2, 3,..., 15) are 2.0744, 4.1469, 
6.9070, 10.3517, 14.4772, 19.2792, 24.7534, 30.8952, 
37.7002, 45.1633, 53.2789, 62.0396, 71.4347, 81.4347 and 
they agree with the numbers given in Ref. |2^. Though 
we have calculated about 18% of eigenenergies we have 
obtained only a small fraction of low- lying states. For 
example, we have not determined energies of the first 
excited state for each 5 < 15. Hence, only a very rough 
approximation of the low-temperature specific heat is cal- 
culated from the first four multiplets (5 < 3). The best 
fit is achieved for (E (l)-E )/k B w 19. 2K, as in Ref.|^. 
Since Ai has been calculated to be about 0.69169 we ob- 
tain J/k B w 27.76 K (the value of 28.8 K given in Ref. || 
is obtained from Lande's rule, i.e. for Ai = |). 

The quotient Eq/J gives the total spin correlation of 
the nearest neighbors. To compare it with other spin 
systems it has to be divided by the maximum eigenvalue 



of a product ns,- • s.+i, i.e. by 52.5 for s 



Hence we obtain Eq /6Js(s - 
quite strong correlation. 



6. 



1) —0.83685, what means 



V. SPIN CORRELATIONS AND THE 
CLASSICAL ORDERING 

The thermodynamic quantities have been already de- 
termined for the Feg cluster and the xesu.lts have been 
compared with the experimental dataSEjI^j Therefore, 
in this paper we focus our attention on the eigenstates 
properties: the spin correlations cu* = (5Zj=i s j s j+r)/® 
for r = 1, 2, 3 and the probability pn of finding the system 
in the Neel state. The latter one is, of course, equal to 
a Ni = a N2> i- e - t° the square of the coefficient of |N1) or 
|N2) in a state under consideration. It has to be stressed 
that oni is nonzero if 5 is even (odd) and T — A\ (B 1; 
respectively). The value of j?n can be used to compare 
quantum states determined with the classical ordering of 
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FIG. 1: Correlations vs energy per spin for S — 4, F = B\. 

antifcrromagnets. 



A. Spin correlations 

In the case 5 = 0, due to isotropy in spin-space, the 
z-correlations of the nearest neighbors, ojf, should be one 
third of energy (per spin) in a given state. This has-been 
confirmed by the calculations presented elsewhere£!3 For 
other states with M = (r = A x , B x , Q < S < 15) 
we have plotted u>* versus energy per spin £ — E/6J. 
Typical drawings are presented in Fig. |l| for S = 4 and 
r = B\. In all cases u>f is almost linear function of £ 
(with a positive slope) and only a few values for high ener- 
gies are larger than zero. On the other hand, wf decreases 
with growing £, but only some low- lying states show pos- 
itive (ferromagnetic) cjf with negative wf and a>|, the 
feature characteristic for the ground state of an anti- 
ferromagnetic system. However, it might happen that 
for u>% < spin correlations of the next-nearest neigh- 
bors are positive for other quantization axes. It could 
be checked considering terms J2j s j ' s i+2 7 what has not 
been done in our calculations. We tried to fit some sim- 
ple functions to but the results obtained occurred 
rather ambiguous. It seems that an exponential function 
wf (£) + !. 132 = exp(-1.224-0.432£) fits better than the 
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FIG. 2: (a) A slope a(S) and (b) an intercept b(S) as a func- 
tion of the total spin S; 'O' and '+': results of linear regression 
ui\(£) = a(S)£ + b(S); lines give (a) linear regressions of ao(S) 
(dashed line), 01 (S) (dotted line) and (b) a parabolic fit to 
b(S) (see text for details). 



power law wf (£) + 3.191 = 6.850 • (£ - £ )- Q - 532 (5 = 4, 
r = Bi). The correlations with the third neighbors (£) 
show positive maximum for all 5 < 14, r = A\ and 
S < 13, r = B\. It is interesting that single states, 
obtained for pairs (S,T) = (13, Bi), (14, B x ), (15, A x ), 
have negative: —0.1833,-0.6327,-0.2155, respec- 
tively; there are no states for pairs (14, A\) and (15, -Bi). 

The correlation uo\ has been analyzed in two aspects. 
At first, assuming linear dependence wf(£) = a(S)£ + 
b(S) for each 5, we have plotted coefficients a(5) and 
6(5) versus 5 (see Fig. ||). Results of linear regressions 
are denoted by the crosses '+' are obtained for pairs 
(5, r) containing the ground state of an 5-multiplet, i.e. 
5 even, T = B\ and S odd, T = Ai\ the diamonds ! 0' 
denote results for the other pairs. Since 0(15*) behaves 
differently in those cases, then the linear functions ao(5) 
(the first case, dashed line) and a x (S) (the other case, 
dotted line) have been considered separately. The points 
for 5 = in both cases and for S = 1 in the first case 
have been omitted due to highly irregular behavior (see 
also below and Fig. ||). As the results we have obtained 

a (5) = 0.02215 + 0.291 , a x (5) = 0.02495 + 0.245 . 

It has to be stressed that ao(5) deviates from a line for 
high 5 and, maybe, a parabolic fit should be applied 
in this region. The intercept 6(5) does not show such 
irregularities so all results are fitted to one (second order) 
function. The best-fit line is given as 

6(5) = -0.01975 2 + 0.05915 - 0.0518 . 

Let u)q(S) denote ojf in the state with the lowest en- 
ergy for each 5-multiplet. The plot in Fig. || shows this 
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FIG. 3: Correlations in the state with the lowest energy in 
each S-multiplet vs the total spin S (diamonds). The dotted 
line: Lande's rule for correlations given by Eq. (0). 
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FIG. 4: The Neel probability pn versus energy £ for states 
with S = 4, T = Bl 



function fitted by a relation analogous to Lande's rule 



u (S) = ±S(S + l) + ±E . 



(6) 



Substituting Lande's rule As = S(S + l)/3 one obtains 



u (S) = l(£o(S) + £ ), 



(7) 



Lande's rule for correlations is broken for S = 0, 1. In the 
first case Eq. @ yields ujo(S) — \£q, whereas the exact 
value is £q/3. The latter case, S — 1, shows deviation in 
the other direction: calculated value of wo(l) is —4.172, 
whereas from Eq. (Q) one obtains —3.632, so the real 
value is about 15% larger in magnitude. Similar results 
have been obtained for six spins \ and |. 



B. Neel probability 

It is interesting to compare the results obtained with 
correlations in the Neel state. In this case z-correlations 
(per bond) are equal to ±6.25 since the Neel state ex- 
hibits the long-range order. The ground state determined 
here has alternating, but decreasing in magnitude, cor- 
relations. However, the total nearest neighbors correla- 
tion, equal to the ground state energy 3wo(0), is larger 
than in the Neel state, since in the Neel state x- and y- 
correlations are equal to zero. The orbit containing the 
both Neel states gives one state with symmetry r = B\ , 
namely |Neel). The absolute value of the coefficient of 
this state in the ground state is = 0.3284, so the 
probability is equal p^ — 0.1078. The probability of a 
single Neel state is pm = pn2 = ^Pn = 0.0539. Such 
small number does not explain so strong nearest neigh- 
bor correlation in the ground state. However, one has 
to take into account that the Neel state, |N1) or |N2), 
represents only one possibility of such spin arrangement 
that each of sublattices (of a bipartite antifcrromagnet) 
has the maximum total spin S max (in the other words — 
each sublattice is ordered ferromagnetically), but mag- 
netizations have opposite directions. For a given spin 



number, S'max in the case considered, there are 2S' max + 1 
possible spin projections, so there are 2S' max + 1 differ- 
ent Neel statesrj each of them corresponds to different 
choice of quantization axis and its orientation. In the 
case considered S'max = 3x1 = 4^, then there are 16 
Neel states. Two of these states are given by |N1) and 
|N2), whereas the others can be obtained as linear com- 
binations of states with S = M — 0. Hence, the global 
probability is p^ = 16p N i = 0.8624. 

To have a bit deeper insight we have analyzed the Neel 
probability p^ for other 5-multiplets. As it has been 
mentioned earlier only states with S even, F = B\ and 
r = A\ for S odd are taken into account. The lowest 
lying state (with energy £q(S)) has always the maximum 
Pn, much higher than in the other states. A rough es- 
timate shows exponential decreasing. For example, in 
the case S — 4, T = B\ the best-fit line is given as 
p N {£) = exp(-2.176£- 14.700) (see also Fig. |). Fitting 
a function exp(a(S')£ + (3(S)) to p^ in each S-multiplet 
we have obtained series a(S), P(S) which can be approx- 
imated by exponential functions a(S) — — exp(0.463S' — 
4.781) - 2.122 and 0(S) = exp(0A63S - 3.061) - 14.643, 
respectively (see Fig. ||). Points for S = 0, 1 have not 
been taking into account due to their irregular behavior. 

All above considerations concern the Neel probability 
Pn- However, as it has been mentioned discussing the 
ground state, the Neel states |N1) and |N2) are only two 
of ns+1 possible states \Sa, m, Sb, —m) with Sa — Sb — 
ins and \m\ < ^ns, where Sa (Sb) denotes the total spin 
in one sublattice of a bipartite antifcrromagnet. Coupling 
of sublattice spins leads to (ns + l) 2 states with the total 
spin < S < ns. In each S-multiplet there is one state 
\S,0) with the magnetization equal to zero. This state 
is a linear combination of ns + 1 states \Sa, fn; Sb, —n%), 
namely 



ns/2 
m—~ns/2 



ins ins S 
m — m 



|ins, m; ins, —m) , 



where the expression in brackets is the Clebsch-Gordan. 
coefficient (CGC) for coupling spin representations _D S .EJ 
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a(S) 



0(S) 




FIG. 5: Coefficients a(S) (squares) and f3(S) (diamonds) of 
estimations pn = exp(a(S)£ + /?(£)) for < S < 13 and the 
best-fit exponential lines (dotted). 



In the case 5 = the appropriate CGC's are equal to 
(_l)i«s-m( ns + l)4 , 



so all states 



^ns, - 



have the same probability 



and, therefore, the global Neel probability p^ = (ns + 
1)pn- Note the Neel states |N1) and |N2) correspond 
to m = ±^ns, respectively. Hence, these states enter 
1 0, 0) with opposite signs for odd ns, according with the 
Marshall rule.C^EJ 

In this paper we consider states with total magneti- 
zation equal to zero and the Neel states correspond to 
sublattice states with = Sb = \m\ = ns/2. So, we 
are interested in the CGC's for.-these values; they can be 
derived from a general formulaEI and one obtains 



«N1(2) 



11c 
±hns =F^ns 



(±1) 



ns+S 



[(ns- S)\(ns + S+ 1)!] 



(8) 



For example, the maximum total spin S = ns yields equal 
coefficients for the both Neel states 



£*N1 = «N2 



(ns)l 



this state always has the symmetry Ai (Tq in a general 
case). To obtain the global Neel probability we have 
to divide ??ni(2) by a^. The results for 5 = 14 and 
5 = 15 can be easily predicted. In both cases there is 
the unique state with these labels and M — 0, so we 
have to obtain |<zni| = \/2|Q!Ni|- In fact, our calculations 
confirm this relation, so in these cases pj, = 1. The global 
Neel probabilities for the states with the lowest energy 
£o(S) in each 5-multiplet are collected in Table III. 

Assuming (as Iske and Caspers in Ref. ^) that the 
global Neel probability is related to the long-range or- 
dering we come to a conclusion that it is present in the 
states \ns,0) and \ns — 1,0). In the other states \S,Q) 



TABLE III: The Neel probability p m , the CGC for the Neel 
state in 15,0), and the global Neel probability p^ = a^ 2 pNi 
for the lowest lying states in each S-multiplet. 



s 


10 3 -p N l[%] 


aT 2 

Nl 


Pn[%] 





5391.126 


16 


86.258 


1 


14305.617 


272/45 


86.470 


2 


18634.880 


816/175 


86.892 


3 


17979.305 


15504/3185 


87.520 


4 


13999.674 


5168/819 


88.340 


5 


9062.557 


15504/1573 


89.324 


6 


4928.709 


15504/845 


90.432 


7 


2254.612 


118864/2925 


91.621 


8 


863.2259 


6992/65 


92.857 


9 


273.7878 


31280/91 


94.111 


10 


70.76512 


594320/441 


95.368 


11 


14.54008 


46512/7 


96.613 


12 


2.286329 


1069776/25 


97.834 


13 


0.2584271 


3447056/9 


98.979 


14 


0.0186955 


5348880 


100.000 


15 


0.00064467 


155117520 


100.000 



this ordering is slightly 'broken' or, in the other words, 
the ground state of quantum antiferromagnetic rings are 
ordered to some extent only. In our case it is reflected 
by the value p^ = 0.86258. However, this number has 
to be compared with p^ dis in a completely disordered 
stateS One of possible definitions of such a state, ap- 
propriate in our considerations, is to combine all vectors 
with M = and T = B\. There are 385 such vectors, so 
a N = 385-5 w 0.051; then p^ dis « 0.021. Much larger 
global Neel probability we obtain taking into account 
only states with 5 = 0. There are twenty such states and 
we take a sum of Neel orbit coefficients divided by \/20 
as «n is this state. In this approach p^ dis w 0.235. The 
value obtained is almost four times smaller than that in 
the ground state. Of course, this ratio should decrease for 
larger rings, since — as it is well known — no long-order ex- 
ists in linear antiferromagnets. On the other hand, how- 
ever, small antiferromagnetic nanoclusters should exhibit 
a well-ordered ground state. 



VI. FINAL REMARKS 

The small ferric wheel Fe6 has been recently inves- 
tigated by many authors, so we have decided to an- 
alyze some less discussed properties — spin correlations 
and probability of finding the system in the Neel state. 
In both cases we obtain results suggesting quite well or- 
dered ground state with highly correlated spins. For ep. 
ample, ujf = -2.441, wf = 1.959, and w| = -1.952.E2 
It should be stressed that z-correlations reach maximum 
magnitude in states with 5=1 (cf. Fig. |J); the values 
are —4.172, 3.907, and —3.897, respectively. However, 
these states are not isotropic, so x- and y-correlations are 
much smaller, whereas total spin correlations Lu r = 3lu* 
for 5 = 0. As regards the Neel probability, deviation 
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FIG. 6: The global Neel probability pjj versus the total spin 
S (diamonds) and the best-fit line (see text for details). 



from the classical ordering in the ground state can be 
observed. On the other hand, = 1 in the cases 
S = 14, 15 shows that these multiplets behave classi- 
cally. It has to be stressed that states with M = 
are considered only. Hence, we obtain energies for all 
states in a given 5-multiplet (for T = A\ , B\ ) . How- 
ever, correlations calculated are related to states \S, 0) — 
they cannot be used for states with M ^ 0. For ex- 
ample, in the case of the unique multiplet S — 15 we 
have uj* = 6.25, r = 1,2,3, for M = ±15, whereas for 
M = we have w z r = -0.2155, r = 1,2,3. Plotting 
pj£j versus the total spin S one can observed almost lin- 
ear dependence for high S (omitting point S = 15). A 
better approximation has been obtained for a function 
p^(S) = exp(-0.31S + 2.84v / S- 8.25) + 0.86 (cf. Eq. (§) 
and Fig. J]). Having determined such relations (see also 
Figs, Q, |^) for ferric dimers Fe2 and the ferric wheel 
Feio we would try to approximate results foiL-Feis. The 
ring Fes recently studied by Waldmann et al.o may lead 
to inconsistent results, since the ground state corresponds 
to r = A% (see Ref. ^ and the end of Sec. III). Note that 
there is lack of Fei4 , but nevertheless it seems interesting 
to investigate such a molecule. There are about 8 • 10 10 
basic states, but 'only' about 5 • 10 9 of them have M = 0; 
the total spin S = leads to 56 267 133 states (there are 
not basic states, however). Taking into account symme- 
try properties will decrease these numbers significantly. 

The methods proposed by Bonner and Fishercd almost 
forty years ago have been mainly used to estimate prop- 
erties of magnetic materials in the thermodynamic limit. 
On the other hand, magnetic nanoclusters have given 
rise to more detailed investigations of small spin systems 
comprising ions with high spin number s. The numeric 
complexity is similar: there are 60466 176 basic states in 
the case of ten spin s = | what is comparable with that 
number for 26 spins s = |. However, our aims are a bit 
different in these cases. Trying to determine properties in 
the thermodynamic limit we observe their dependencies 
on number of spins n, whereas investigating nanomagnets 
we rather try to determine a 'best-fit' model of a given 
molecule. Therefore, in the first case the considerations 



are mainly limited to one, sometimes simplified, model, 
while fitting to experimental data we may need to inves- 
tigate some different and rather complicated models. 

Amongst many magnetic molecules those containing 
Fe 3+ ions seem to be very interesting, since theji-.can be 
used as models of such metalloprotein as ferritin£j A fact 
that this molecule contains as many as 4500 Fe ions does 
not depreciate validity of results obtained for molecules 
comprising several magnetic ions. Experimental data 
show that noncompensate spin in ferritin pis about 50 
and it originates from surface randomness .EJ Moreover, 
these data suggest that ferritin exhibits quantum tunnel- 
ing similar to that observed in Mni2 molecules. 

Since recently some, pap ers using a similar approach 
have been publishecOoE3 we would like to stress differ- 
ences and advantages of the method applied in this paper. 
At first, we go much deeper into structure of the permu- 
tation representations introducing additional indices (or- 
dered and nonordered partitions) for basic states. This 
also leads to some simplifications in the procedure gen- 
erating operatop-patricesJ^I'Lj A general formula for ma- 
trix elementsHEa significantly reduces a number of cases 
to be considered. Note that this formula needs introduc- 
tion of double cosets, a structure not mentioned in the 
earlier papers. Moreover, many results for the dihcdrjal 
groups can be presented in easy-to-use analytical form.Ej 
We do not need to use any projection procedure to de- 
rive irreducible basis, because it is determined for tran- 
sitive, representations of the dihedral groups in a general 
form.EjCj Using » package for multiple precision arith- 
metic, like GMP,EJ we are able to find exact eigenvectors 
of S 2 operator in the case of one-dimensional irreps (and 
also for two-dimensional irreps in some simple cases, e.g. 
for Dg or Ds, when cos(2n/n) can be expressed by square 
roots of integers) . This allows determination of spin cor- 
relations and Neel probabilities with high accuracy. 
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APPENDIX A: AN EXAMPLE 



To illustrate the method for generating operator ma- 
trices the case of four spin s = I is presented below. The 
solutions of eigenproblems of S 2 and Ti. as well as dis- 
cussion of the spin correlations and Neel probability are 
left out, since at the actual stage of the project they are 
performed by means of (more or less) standard numeric 
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methods (the only change is application of procedures 
from the GMP library). 

Assuming the isotropic Heisenberg Hamiltonian both S 
and M are good quantum numbers. There are 4 4 = 256 
Ising configurations and the maximum total spin num- 
ber is S = 6; hence, the possible magnetizations are 
M = 0,±1,...,±6. The dimensions of subspaces Lu 
containing states with a given M are as follows: 



M 


±6 


±5 


±4 


±3 


±2 


±1 





dim L m 


1 


4 


10 


20 


31 


40 


44 



It allows determination of numbers n(S) of S-multiplets 



for each S as n(S)f^ dim L 



jM=s ' 



dim L 



and ra(S = 6) = 1;E3 hence one obtains: 



M=S+1 



for S < 6 



s 


6 5 


4 


3 


2 1 


n(S) 


1 3 


6 


10 


11 9 4 



Of course we have ^]j =0 (2S + l)n(S) — 256 and 

Es=o n (^) = 44 = dimLjv/=o since each S-multiplet 
contains the state \S M = 0). To investigate nonmag- 
netic states with M = one has to decompose the 44- 
dimensional space Lq into seven subspaces labeled by S 
and with dimensions given in the table presented above. 

To begin with we start from considering ordered par- 
titions of n = 4 into no more than four parts. There are 
five such partitions, namely 

[4], [3,1], [2,2], [2,1,1], [1,1,1,1]. 

Each of this partition represents a class of nonordered 
partitions of n = 4 into four non-negative parts ko,ki,k2, 
and fca. There are 35 such partitions collected in Ta- 
ble IV. From this table one can learn that there are five 



nonordered partitions with y\ r 



Ik, 



6, i.e. the mag- 



netization M = (see Eq. (0)). These five partition 
come in three types: [3,1], [2,2], and [1,1,1,1]. There- 
fore, to investigate states with M — one has to con- 
sider only these three ordered partitions. The number 
of configurations of type [k] is determined by the poly- 
nomial coefficient n!//co!fci!fc2!&3! and is the same for all 
nonordered partitions of a given type [«;]. In the case 
considered, i.e. for [«] = [3, 1], [2,2], and [1, 1, 1, 1] these 
numbers are 4, 6, and 24, respectively, then there are 
four states corresponding to [k] — [0, 3, 0, 1], [1, 0, 3, 0], 
six states for [k] = [2, 0, 0, 2], [0, 2, 2, 0], and 24 states for 
[k] = [1, 1, 1, 1]; 44 states in total. For example, the par- 
tition [1, 0, 3, 0] labels a four-element set of states: 



{|- 



since relations ko = 1, &2 = 3, and k\ 



means 



that there are one spin projection — | and three projec- 
tions \. On the other hand, the unique orbit labeled by 
[k] = [1,1,1,1] contains 4! = 24 states with all projec- 
tions different: |-§, -|, |, §), |-| , -f, |, |) etc. 

The decomposition of an orbit 0[k] (of the symmetric 
group £4) into orbits of a Hamiltonian symmetry group 



TABLE IV: Classes of nonordered partitions [k] labeled by 
ordered partitions [k] for n = 4 and s = §; in the third 
column the magnetization M = 5Zf = o ^1 ~ 6 is calculated. 



[«] 


[k] 


M 


orbit representative 


[4] 


[4,0,0,0] 


-6 


-3/2,-3/2,-3/2,-3/2) 




[0,4,0,0] 


-2 


-1/2,-1/2,-1/2,-1/2) 




[0,0,4,0] 


2 


+ 1/2, +1/2, +1/2, +1/2) 




[0,0,0,4] 


6 


+ 3/2, +3/2, +3/2, +3/2) 


[3,11 

I ' J 


[3,1,0,0] 


-5 


- 3/2, -3/2, -3/2, -1/2) 




[3,0,1,0] 


-4 


-3/2, -3/2, -3/2, +1/2) 




[3,0,0,1] 


-3 


-3/2, -3/2, -3/2, +3/2) 




[1,3,0,0] 


-3 


-3/2,-1/2,-1/2,-1/2) 




[0,3,1,0] 


-1 


-1/2, -1/2, -1/2, +1/2) 




[0,3,0,1] 





-1/2, -1/2, -1/2, +3/2) 




[1,0,3,0] 





-3/2, +1/2, +1/2, +1/2) 




[0,1,3,0] 


1 


- 1/2, +1/2, +1/2, +1/2) 




[0,0,3,1] 


3 


+ 1/2, +1/2, +1/2, +3/2) 




[1,0,0,3] 


3 


-3/2, +3/2, +3/2, +3/2) 




[0,1,0,3] 


4 


- 1/2, +3/2, +3/2, +3/2) 




[0,0,1,3] 


5 


+ 1/2, +3/2, +3/2, +3/2) 


[2,2] 


[2,2,0,0] 


-4 


-3/2,-3/2,-1/2,-1/2) 




[2,0,2,0] 


-2 


-3/2, -3/2, +1/2, +1/2) 




[2,0,0,2] 





-3/2, -3/2, +3/2, +3/2) 




[0,2,2,0] 





- 1/2, -1/2, +1/2, +1/2) 




[0,2,0,2] 


2 


- 1/2, -1/2, +3/2, +3/2) 




[0,0,2,2] 


4 


+ 1/2, +1/2, +3/2, +3/2) 


[2,1,1] 


[2,1,1,0] 


-3 


-3/2, -3/2, -1/2, +1/2) 




[2,1,0,1] 


-2 


-3/2, -3/2, -1/2, +3/2) 




[2,0,1,1] 


-1 


-3/2, -3/2, +1/2, +3/2) 




[1,2,1,0] 


-2 


-3/2, -1/2, -1/2, +1/2) 




[1,2,0,1] 


-1 


-3/2, -1/2, -1/2, +3/2) 




[0,2,1,1] 


1 


-1/2, -1/2, +1/2, +3/2) 




[1,1,2,0] 


-1 


-3/2, -1/2, +1/2, +1/2) 




[1,0,2,1] 


1 


-3/2, +1/2, +1/2, +3/2) 




[0,1,2,1] 


2 


-1/2, +1/2, +1/2, +3/2) 




[1,1,0,2] 


1 


-3/2, -1/2, +3/2, +3/2) 




[1,0,1,2] 


2 


-3/2, +1/2, +3/2, +3/2) 




[0,1,1,2] 


3 


-1/2, +1/2, +3/2, +3/2) 


[1,1,1,1] 


[1,1,1,1] 





-3/2, -1/2, +1/2, +3/2) 



G (isomorphic to the dihedral group D4 in the case con- 
sidered) depends on its type, i.e. on the ordered partition 
[k]. Before presenting these decompositions a brief out- 
look of the group D4 is necessary. This group is generated 
by the four-fold rotation C4 and the two-fold rotation Uq 
with the generating relations 



Ct = E, U* = E, {C A Uv) 2 = Ul 



E. 



where E denotes the unit element in D4. The labeling 
scheme of group elements and numbering of nodes (spins) 
is presented in Fig. 0. There are five classes of conjugated 
elements: {E}, {C 2 }, {C 4 , C7 1 }, {Uq, U 2 }, {U u U 3 } and, 
hence, five irreducible representations: A\, A%, Bi, B2, 
Ei with characters presented in Table [y| The dihedral 
group D4 has ten subgroups collected in eight classes 
of conjugated subgroups (in this very simple case only 
two classes contain two subgroups). However, three sub- 
groups are excluded from a set of possible configuration 
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FIG. 7: Labeling scheme of elements in the dihedral group D4 
and nodes (spins) for a system of four spins. A black square 
denotes the four-fold axis and thick lines denote two-fold axes. 



TABLE V: Characters of irreducible representations of the 
dihedral group D6. 



L 


{E} 


{C , 4 ,C 4 - 1 } 


{C2} 


{U ,U 2 } 


{Ui,U 3 } 


Ai 


1 


1 


1 


1 


1 


A 2 


1 


1 


1 


-1 


-1 


Bi 


1 


-1 


1 


-1 


1 


B 2 


1 


-1 


1 


1 


-1 


£1 


2 





-2 









stabilizers.E3H It results in five classes of subgroups or, 
in the other words, five types of orbits. These classes are 
represented by subgroups: 

d = {£}, D° = {E,U }, V{ = {E,Ui}, 

D° = {E,C 2 ,U ,U 2 }, 

D 4 = {E, C 4 , C 2 , C7/ 1 , U , Ui, U 2 , U 3 } . 

The four-element orbit labeled by the partition [3, 1] 
is also an orbit of the dihedral group D4 with a stabi- 
lizer Dj. As representatives of orbits of this type one 
can choose a state in the form \a,b,b,b). In the case 
considered here, two orbits of this type are represented 
by states ||, — |, — |, — |) and | — |, |, |, |), respectively. 
The orbits of the type [2, 2] have been considered in the 
main text. Each of them decomposes into two orbits of 
the group D4: the first has a stabilizer D}, whereas the 
second one has a stabilizer D§. In this way one obtains 
four new orbit representatives: 

[2 21 ' I- 3 - -2 2 £\ 13 3 _3 3\. 

1 2 , 2' 2» 2/ ' I 2'2' 2'2/' 

fO 2 2 01 • l-i -i i i) l-i - -- -) 

[u, i, «, uj . 1 2 , 2 ! 2» 2/ ' I 2'2' 2'2/' 

The last, 24-element, orbit decomposes into three orbits 
in the same type: all of them have trivial stabilizer Qi-, 
so each of them is the so-called regular orbit of D4.E3 
As representatives of these eight-element orbits one can 
choose states 

|_3 _1 1 3\ |_3 1 3 _1\ I 3 3 _1 1\ 

I 2' 2' 2 ' %i ' I 2 ' 2' 2 ' 2/ ' I 2' 2 ' 2' 2 / ' 

In this way the complete group-theoretical (and com- 
binatorial) classification of all 44 Ising configurations 



with M — has been done. For example the state 
/i = ||)— 5) is labeled by the nonordered parti- 
tion [k] — [0, 2, 2, 0] (so the ordered partition is [2, 2]); its 
stabilizer is D 2 so it belongs to the orbit represented by 
v = \ — o> 2> — 3' 2) an< ^ as a representative g r of the left 
coset g r D 2 € D4/D2 one can choose the four-fold rotation 
C4, since [i — C/^v. 

The partition [fc] may be used to simplify determina- 
tion of matrix elements of S 2 or H. Assuming that only 
bilinear terms are present one can notice that states la- 
beled by [2, 0, 0, 2] can be transformed into states labeled 
by [1,1,1,1] only. States of the first type comprise two 
projections I and two projections — |. Therefore, the 
only nontrivial action of a bilinear term s^s^ is possible 
when Sj~ acts on the projection — | and acts on |. It 
results in the state containing all four different projec- 
tions, i.e. of the type [1, 1, 1, 1]. Similar discussion of all 
five partitions leads to the following graph:Eilc3 



[0,3,0,1] 




[1,0,3,0] 

The lines join partitions, which are related by the action 
of bilinear terms sfsZ: encircled partitions denote loops 
in this graph, i.e. for each state labeled by these parti- 
tions exist such a term s^s^ that a state labeled by the 
same partition is obtained after its action. For example, 
S + S -|_I I I _ I\ - All _I I _I\ 

°1 °2 \ 2'2'2' 21 12' 2'2' 2'' 

The obtained orbits have four different stabilizers: Ci, 
T>i, D{, and D§. Note that there are three orbits of 
the first type and two orbits of each of the other types. 
To introduce the linear structure one has to decompose 
(transitive) permutatk>p,representation BP i:U into (lin- 
ear) irreducible ones.at 3 ] For these four groups these de- 
compositions are as follows: 

i?° 4:Cl = Ax © A 2 © B 1 © B 2 © 2Ex , 

E D 4 :D? = Ax®Bx®Ex, 
R B 4 :Bl = A 1 ®B 2 ®E ll 

Therefore, the 44-dimensional space Lq can be decom- 
posed into nine subspaces labeled by T = A\, three by 
r = A 2 , seven by T = B\, five by r = B 2 , and ten by 
r = E\. Analogous considerations for 40-dimensional 
space Lm=i give us numbers 7, 3, 7, 3, and 10, respec- 
tively. It means that amongst four (44 — 40) multiplcts 
with S — two (9 — 7) are labeled by Ai, and the other 
two (5 — 3) by B 2 . In the case of antiferromagnetic in- 
teractions the ground state has the symmetry A±, so it 
belongs to the first pair. 
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